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Abstract. Photoproduction of neutral pions from nuclei (carbon, calcium, niobium, lead) has been studied 
for incident photon energies from 200 MeV to 800 MeV with the TAPS detector using the Glasgow photon 
tagging spectrometer at the Mainz MAMI accelerator. Data were obtained for the inclusive photoproduc- 
tion of neutral pions and the partial channels of quasifree single tt", double tt", and tt^tt photoproduction. 
They have been analyzed in terms of the in-medium behavior of nucleon resonances and the pion - nucleus 
interaction. They are compared to earlier measurements from the deuteron and to the predictions of a 
Boltzmann-Uehling-Uhlenbeck (BUU) transport model for photon induced pion production from nuclei. 

PACS. 13.60.Le meson production - 25.20.Lj photoproduction reactions 



1 Introduction 

The study of meson photoproduction from nuclei is moti- 
vated by two strongly interrelated aspects, namely possi- 
ble in-medium modifications of hadrons and the meson - 
nucleus interaction. The in-medium properties of mesons 
and nucleon resonances is a topic which is hotly debated. 
However, with the exception of the Z\-isobar, experimental 
results are still scarce and partly contradictory. In-medium 
modifications of nucleon resonances may arise from many 
different effects. The most trivial one is the broadening of 
the excitation functions due to Fermi motion of the bound 
nucleons. Pauli-blocking of final states decreases the decay 
widths but this is counterbalanced by coUisional broaden- 
ing due to resonance-nucleon scattering. In particular for 
the Di3(1520) it has been found in T,'5"3' that a large 
collisional broadening can be generated from tt and p ex- 
change and that the in-medium width is sensitive to modi- 
fications of the spectral functions of these mesons. Besides 
such conventional effects arising from meson-nucleon scat- 
tering, additional mechanisms [5||2 could lead to a shift 
of spectral strength of the p meson to smaller invariant 
masses, thus further enhancing the width of the D13. 

Photon induced reactions are particularly well suited 
to study such effects since photons probe the entire nu- 
clear volume. Reactions induced by strongly interacting 
particles like pions are always restricted to the surface re- 
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gion and thus probe a lower effective density. The first 
experimental investigation of the nuclear response in the 
resonance region to photons was done with total photoab- 
sorption measurements jHIIZl for nuclei ranging from ^Li 
to ^•^^U. Total photoabsorption has the advantage that no 
final state interaction effects (FSI) must be accounted for 
so that indeed the entire nuclear volume is probed. As a 
consequence, when scaled by the atomic mass number, the 
total cross sections from all nuclei fall practically on the 
same curve ('universal curve'). This curve shows a clear 
signal for the A resonance but no evidence of the bumps 
observed for the free nucleon in the second and third res- 
onance regions (see fig.DJ. This result of the experiments 
was surprising and is even today not fully understood. 

The situation is complicated by the already complex 
structure of the resonance peak for the free proton. The 
resonance consists of a superposition of reaction channels 
which differ in their energy dependence (see fig. 01. Most 
of the rise of the cross section to the maximum around 700 
MeV comes from double pion production. The interpreta- 
tion of double pion photoproduction is complicated (see 
_8 for a recent overview), but it is known that contribu- 
tions, which are not related to the excitation of resonances 
in the second resonance region, are significant. In particu- 
lar, Z\-Kroll-Rudermann and A pion-pole terms dominate 
the TT+TT^ final state which has the largest cross section 
of all double pion channels. Gomez Tejedor and Oset [H] 
have pointed out that for this channel the peaking of the 
cross section around 700 MeV is related to an interference 
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between the leading Z\-Kroll-Ruderniann term and the se- 
quential decay of the D13 resonance via D13 -^ Att. Hirata 
et al. ^ni have subsequently argued that the change of 
this interference effect in the nuclear medium is one of the 
most important reasons for the suppression of the bump 
structure. 
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Fig. 1. Total photoabsorption on the proton and on nuclei. 
The symbols show the total cross section and the decomposi- 
tion into partial reaction channels for the proton. |llll2lli-{ll4l 
I15I16| . Solid curve: average of the total photoabsorption cross 
sections for nuclei scaled to their mass numbers [Tj. 



Inclusive reactions like total photoabsorption alone do 
not allow a detailed investigation of such effects. A study 
of the partial reaction channels is desirable. The experi- 
mental identification of exclusive final states is more in- 
volved and FSI effects must be accounted for. The in- 
terpretation of exclusive measurements therefore always 
needs models which account for the trivial in-medium and 
FSI effects like absorption of mesons and propagation of 
mesons and resonances through nuclear matter. On the 
other hand, as a by-product, the analysis of the FSI ef- 
fects enables the study of meson-nucleus interactions. 

During the last few years, in a series of experiments 
with the TAPS detector at MAMI in Mainz, we have 
studied exclusive photoproduction reactions from heavy 
nuclei which can be related to the excitation of certain 
resonances or production mechanisms. The most specific 
reaction in this respect is the photoproduction of 77- mesons 
which proceeds in the second resonance region almost en- 
tirely through the excitation of the Sii(1535) resonance 
| 13II17| . In this case no unexplained depletion of the in- 
medium strength was found |18| , which is in line with theo- 
retical findings .3-19, that the change of the Sn self energy 
in the medium is rather small. The data were in excellent 
agreement with model calculations that take the trivial 



in- medium effects and FSI into account |2()ll21j . A similar 
result was later found by Yorita et al. P? who measured 
the 77 excitation function to somewhat higher incident pho- 
ton energies. In a recent study [221 j i* '^^^ shown that the 
data could be described over the full energy range by ap- 
plying a momentum dependent Sn nuclear potential. An 
attempt to study the in- medium properties of the D13 res- 
onance was undertaken with a measurement of quasifree 
single tt" photoproduction [22 ■ Again in contrast to total 
photoabsorption the second resonance bump was clearly 
visible. On the free nucleon, the structure in this reaction 
is almost entirely due to the D13 resonance. A compari- 
son to deuteron data showed no broadening in nuclei and 
provided no evidence for an increased in-medium width 
of the D13. On the other hand, theoretical models [312] 
predict such a modification. However, as discussed below, 
exclusive reaction channels are dominated by the nuclear 
surface region where in-medium effects are smaller. This 
is mostly due to FSI effects. Furthermore, as discussed 
in [23, even in the absence of FSI, resonance broaden- 
ing effects are suppressed in reactions which do not con- 
tribute to the broadening. The reason is, that when the 
total width r' of a resonance increases as function of the 
nuclear density p (e.g. due to coUisional broadening), the 
branching ratio bf — Ff/F for any other decay channel /, 
which does not contribute to the broadening, decreases as 
function of density. Due to this effect, contributions from 
the broadened resonance to such exclusive final states are 
suppressed in the average over the nuclear volume. 

In the present paper we collect the results for all reac- 
tion channels which involve 7r° photoproduction from nu- 
clei. Data have been obtained for the inclusive reaction, 
where inclusive means that all final states with at least one 
7r° meson were accepted. Furthermore, the quasifree con- 
tributions of single 7r°, double 7r°, and tt^tt^ photoproduc- 
tion have been determined. The data are compared to the 
results of a semi-classical Boltzmann-Uehling-Uhlenbeck 
(BUU) transport model [^^HT^ which accounts for medium 
effects such as Fermi motion, nuclear binding, Pauli block- 
ing, resonance broadening and includes a realistic descrip- 
tion of the FSI. The treatment follows the trajectories of 
the particles through the nucleus and therefore yields ad- 
ditional information not accessible to experiment. 



2 Experimental setup 

The experiments were carried out at the Glasgow tagged 
photon facility 27 installed at the Mainz MAMI acceler- 
ator [^ . The tagged photon facility uses bremsstrahlung 
photons produced with the 855 MeV electron beam in a 
radiator foil. The scattered electrons are detected in the 
focal plane of the tagging spectrometer with an energy res- 
olution of «2 MeV. The experiments covered the photon 
energy range from 200 - 800 MeV. Solid targets with thick- 
nesses of 25 mm (C), 10 mm (Ca), 2.1 mm (Nb), and 0.47 
mm (Pb), corresponding to roughly 0.1 radiation lengths 
were irradiated. The diameter of the targets was 50 mm 
and the beam spot size at target position was approxi- 
mately 30 mm. The neutral pions were detected, via their 
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in detail in |16| for the measurement with the deuterium 
target. 



Fig. 2. Setup of the TAPS detector at the Mainz MAMI accel- 
erator. The beam entered the target chamber from the lower 
right edge. 



two photon decays, with the electromagnetic calorimeter 
TAPS [58^,'2cr. The individual modules of the calorime- 
ter are hexagonally shaped BaF2 crystals of 25 cm length 
with an inner diameter of 5.9 cm. The inclusive and sin- 
gle tt" photoproduction data were obtained with the setup 
discussed in 16 . The BaF2 modules were arranged in 5 
blocks, each containing 8x8 modules. With the exception 
of the backward block, all modules were equipped with 
individual plastic veto detectors. The blocks were placed 
in one plane around the target at a distance of 55 cm at 
polar angles of ±38°, ±88°, and +133°. The double pion 
production reactions were studied with an upgraded ver- 
sion of the TAPS detector covering a larger solid angle 
(see fig. 12) • In this setup 6 blocks with 64 modules and 
a forward wall with 138 BaF2 modules were used. The 
blocks were placed 55 cm away from the targets at polar 
angles of ±54°, ±103°, and ±153°, while the forward wall 
was placed 60 cm away from the target center. 



3 Data analysis 

The discrimination between photons and particles was 
done with the aid of the plastic veto detectors, a time- 
of-flight analysis, and a pulse shape analysis of the signals 
from the BaF2 scintillators in the same way as discussed in 
|16) . Neutral pions were identified with a standard invari- 
ant mass analysis of photon pairs. In addition to the anal- 
ysis in 221 , charged pions for the extraction of the tt^tt^ 
cross section were identified from time-of-flight versus en- 
ergy plots. A typical example for this analysis is shown 
in figure 13 Photons have already been eliminated by the 
methods mentioned above and only charged hits, with an 
identified tt" meson in coincidence, are included. Back- 
ground from random coincidences between TAPS and the 
tagging spectrometer was determined by shifting the time- 
of-flight cut away from the prompt peak. The absolute 
normalization of the cross sections was done as described 
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Fig. 3. Time-of-flight versus the energy deposited in the BaF2 
detectors. Three bands corresponding to deuterons, protons, 
and charged pions are visible (photons already eliminated). 



3.1 Inclusive n° photoproduction 

All events with at least one n'^ identified via the invari- 
ant mass analysis were accepted without any further kine- 
matic cuts. Combinatorial background from pairs of pho- 
tons originating from the decay of different pions was fit- 
ted and subtracted in the invariant mass spectra (see ref. 
[TB] for typical invariant mass spectra). For differential 
cross sections this was done separately for each bin of the 
respective observable. 

The energy and angle dependent detection efficiency 
of the TAPS detector was determined with Monte Carlo 
simulations, carried out with the GEANT3 code jSO]. It 
was previously shown that such simulations reproduce the 
photon response of TAPS precisely J5S]. Inclusive 7r° pro- 
duction is a complicated mixture of different reaction chan- 
nels. With the exception of coherent single tt'^ production 
all channels have at least a three-body final state and 
all are subject to final state interaction processes. Con- 
sequently, no reliable models exist which could be used as 
event generators for the Monte Carlo simulations. How- 
ever, the detection efficiency e{0^a, T^o) can only depend 
on the laboratory polar angle 0^o and the laboratory ki- 
netic energy Tt^o of the pions (the dependence on the az- 
imuthal angle <P can be easily integrated out since no po- 
larization degrees of freedom are involved). Since O^^o and 
T^o are measured, the efficiency can be corrected indi- 
vidually for each detected 7r° meson. This method gener- 
ates a model independent result for the cross section as 
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Fig. 4. Detector acceptance for vr -mesons as a function of 
the laboratory polar angle and the laboratory kinetic energy 
of the mesons. Upper part: efficiency calculated with a Monte 
Carlo simulation. Lower part: measured distribution of pions, 
not corrected for acceptance, for incident photon energies from 
200 - 800 MeV. 



long as the efficiency does not vanish for any kinemati- 
cally possible combination of O^o and T^^a . The simulated 
efficiency is plotted in the upper part of fig. ^ The plot 
demonstrates, that although the overall efficiency is not 
high, there is only a small region in the {Ot^o^T^o) plane 
at kinetic energies larger than 200 MeV and forward an- 
gles smaller than 15° were no 7r° mesons can be detected. 
This is so, since in the energy range of interest almost al- 
ways some fraction of the decay photon pairs of t:^ mesons 
falls within the spectrometer acceptance. Kinetic energies 
larger than 700 MeV are impossible for incident photon 



energies <800 MeV and the zero-efficiency region at back- 
ward angles also corresponds to a kinematically forbidden 
region. As an example of a measured distribution the bot- 
tom part of fig. 01 shows the data for ^°Ca. The small ac- 
ceptance gap at forward angles was corrected with a phase 
space simulation. This is not model independent, but the 
largest corrections to the total cross sections are below 5% 
so that possible systematic uncertainties arising from the 
acceptance gap are certainly below the 1% level. 

The inclusive cross section Uinci determined this way 
can be written in the following form: 

<^mcl — O-^f + 2o-27rO + (^■k°tt+ + (^tt°-k- + P<^ri (1) 

(3 — 36,,_»3^o + 6^^^o^+^- « 1.2 

where ax are the partial cross sections labeled with the 
different final states and bi are the branching ratios of the 
77 meson for the respective decay channels. It was previ- 
ously shown in jlB) that this relation holds very accurately 
in case of the deuteron, where all partial cross sections are 
known. This means in particular that final states with pion 
multiplicity larger than two, which do not originate from 
decays of the rj meson, are negligible in this energy range. 
Provided the partial cross sections for 77 and double 7r° 
production are known, the part anm of the total photoab- 
sorption cross section aabs with at least one neutral pion 
or one 77 meson in the final state follows from: 



cl - ^2^0 - (/3 - l)cr^ 



(2) 



The total photoabsorption cross section aabs can then be 
split into: 

O'abs — <^nm + O'ctn + <^r (<jj 

where acm corresponds to reactions which have no neutral 
meson but at least one charged pion in the final state, and 
cr,. corresponds to final states without any mesons. In case 
of the deuteron, ar equals the deuteron photo-breakup. 



3.2 Exclusive reaction channels 

The decomposition of the inclusive process into exclu- 
sive reaction channels is discussed in detail in |lfij for the 
deuteron target, but it is less straight forward for heavy 
nuclei. The larger Fermi momenta cause a stronger smear- 
ing of kinematical observables, and FSI affects the reac- 
tion products. The first step is the separation of quasifree 
single tt" photoproduction. Since the solid angle coverage 
of the detector amounts only to «25% of Air this must be 
done with cuts on the reaction kinematics. Missing energy 
(AEj^) spectra can be constructed for the kinetic energy 
E* of the pion in the center-of-momentum (cm) frame of 
the incident photon and a stationary nucleon: 



AE^ ^ EUja2) ~ EtiEb) 



(4) 



Here, E*(Eb) is the pion energy derived from the energy 
of the incident photon beam Eb and E*{-f 1J2) is recon- 
structed from the momenta of the pion decay photons. 
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Fig. 5. Missing energy spectra for vr" photoproduction for in- 
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cident photon energies 395 - 440 MeV calculated under the Fig. 6. Missing energy spectra for tt" photoproduction for in- 
assumption of quasifree kinematics. cident photon energies of 790 MeV calculated under the as- 

sumption of quasifree kinematics. 



The reaction kinematic is assumed to be quasifree, i.e. 
the pion is produced from one participant nucleon and the 
rest of the nucleus acts as spectator. The neglect of the 
Fermi momenta of the bound nucleons simply broadens 
the missing energy peaks. At low incident photon energies 
a significant fraction of tt*^ photoproduction stems from 
the coherent reaction, in which the recoil is not taken by 
a nucleon but by the whole nucleus. A previously pub- 
lished analysis |3I] has separated the coherent parts with 
a missing energy analysis assuming coherent kinematics 
which gives much narrower peaks since no Fermi smearing 
is involved. In the present analysis contributions from the 
coherent reaction fall within the broader peak of quasifree 
pion production. 

An example of the missing energy spectra at medium 
incident photon energy around 420 MeV is shown in fig. 
There is a striking difference between the missing en- 
ergy distribution from the deuteron and from heavy nu- 
clei. At these energies, in case of the deuteron a peak with 
a small tail at negative values is visible. The tail is due 
to the onset of double pion production which has, how- 
ever, a relatively small cross section in this energy region. 
On the other hand, the distribution for lead shows a sub- 
stantial non-quasifree background. At this photon energy 
the background arises mainly from single pion produc- 
tion with significant final state interaction of the pion or 
from two-body contributions to the production process. It 
must be emphasized that the single n'^ production cross 



section a'' a obtained from this kinematical analysis is not 
identical with the cross section cTttO in eq. ^ By defini- 
tion a^o is the cross section for events with one and only 
one neutral pion in the final state, it does not matter if 
the quasifree energy relation is valid. However, single ir'^ 
events, in which the tt^ suffers an FSI or was produced by a 
two-body interaction, will in general be excluded from ct^q . 
Only for the deuteron, where FSI and two-body absorp- 
tion do not significantly contribute, we have a'^i ~ a^^o. 
At higher incident photon energies, much larger contri- 
butions from multiple pion final states are visible in the 
missing energy distributions. An example is shown in fig. 
|S1 This figure shows also the lineshape of the missing en- 
ergy peaks, which are broader for the heavy nuclei due to 
Fermi motion. The shape of the background (solid curves 
in fig. El was fitted to the deuteron data and then scaled 
for the nuclei. The similarity of the background shapes for 
the deuteron and the heavy nuclei is not surprising. The 
background for the deuteron and the nuclei is dominated 
by the same source, namely double pion events where only 
one pion was detected. Folding the curve with a typical 
momentum distribution for nucleons bound in heavy nu- 
clei has not much effect since already the slope for the 
deuteron is not steep. However, fitting and subtracting 
the background contributions would introduce systematic 
uncertainties which are difficult to control for the com- 
parison to model predictions. We have therefore extracted 
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the quasifree single 7r° production with a cut in the miss- 
ing energy spectra from - 200 MeV, assuming symmetric 
line shapes. Contributions from the tail of the background 
distribution were assumed to be negligible as suggested by 
fig. El The same analysis was applied to the results of the 
BUU calculations so that systematic effects in the com- 
parison between data and model results are minimized. 

The analysis of double tt*^ photoproduction from the 
proton and the deuteron is discussed in detail in |16II15[ 
ESISni- At incident photon energies below the 77 produc- 
tion threshold the cross sections of all other reactions with 
at least three photons in the final state are negligible com- 
pared to double 7r° production. Therefore, an analysis of 
all events with at least three identified photons is suffi- 
cient for the double tt'^ channel. Background may only 
arise from particles which are misidentified as photons or 
from false photon pairs which result from the split-off of 
an electromagnetic shower into two geometrically unre- 
lated components. Particles were completely suppressed 
with the veto detectors, the pulse shape analysis and the 
time-of-flight analysis. Subsequently, a standard invariant 
mass analysis was used to assign two photons (labeled 71 , 
72) to the decay of one tt". The third photon (73) was 
then a candidate for a second 7r° decay with one unde- 
tected photon. Among the latter, split-off photons were 
suppressed with cuts on the relative timing {At-y-y <1 ns 
for all photon combinations), the minimum energy of the 
third photon {E-y.^ >45 MeV), and the minimum opening 
angle between the third photon and the other two photons 

1^7173' ^7273 ^ ^U j. 

The above analysis was only used for the data mea- 
sured with the first setup where the small solid angle cov- 
erage resulted in a very low detection efficiency for events 
with four photons. A much more restrictive analysis is 
used for events with four detected photons which can be 
combined into two pairs with the invariant mass of the 
pion. A two-dimensional representation of the invariant 
masses of the photon pairs is shown in fig. All possible 
combinations of the four photons into two pions are tested 
and the 'best' combination is chosen via a x^ analysis. 
Only these 'best' combinations are shown in the figure. 
After a cut on the invariant mass peak a missing mass 
analysis as discussed in ^^1 serves as a kinematical con- 
straint. The mass of the missing particle (the participant 
nucleon) is calculated from the momenta P~y^ and the en- 
ergies Ej^ of the pion decay photons and the beam energy 
Eb and compared to the nucleon mass rriN. Fermi motion 
is again neglected. The missing mass AM is given by: 



AM 



\ 



{Eb 



rriN 






r-iPt,-j2p^ 



rriM. 



(5) 

Missing mass spectra for calcium below and above the 
77 production threshold are shown in fig. |S1 At incident 
photon energies below the rj production threshold only a 
broad peak around zero is visible which agrees with the 
result of a Monte Carlo simulation taking into account 
Fermi motion, a simple approximation for the final state 
interaction, and the instrumental acceptance and detec- 
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Fig. 7. Invariant mass of photon pair 1 versus invariant mass 
of photon pair 2 for events with four detected photons. 



tion efficiency. At higher incident photon energies back- 
ground from the r] — > 37r° decay appears. In contrast to 
the deuteron case 16,32 , the broadening of the peak due 
to Fermi motion and FSI is so large that no clear separa- 
tion from the rj background is possible. However, since the 
cross section for 77 photoproduction from nuclei is known 
from the analysis of the 27 decay channel of the r] meson 
it can be subtracted. Simulations of double tt" production 
and the 77 background are compared to the data in fig. 
IHl The shape of the sum of the two contributions again 
agrees with the data. 

An analogous analysis was done for the tt^tt^ reaction. 
The relative contribution of the 77 background is much 
smaller in this case. For the free nucleon cr^'-fN — > tt^tt^) 
is almost six times larger than a{'yN -^ tt'^tt'^), while the 
branching ratio for 7/ -^ tt+tt^tt^ is a factor of 1.4 smaller 
than for 77 — > Btt^. 

The detection efficiency for double tt'^ events cannot be 
constructed in the model independent way used for inclu- 
sive and single tt^ photoproduction since the angles and 
energies of the pions are correlated. The case of tt^tt^ pho- 
toproduction is even more difficult since certain regions of 
the polar angle of the charged pions are not within the ac- 
ceptance of the detector (spaces between blocks). There- 
fore a phase space simulation of the reactions, including 
the effects of Fermi motion and a simple approximation of 
meson final state interaction was used. For the latter only 
re-scattering of the pions was considered and the strength 
of this FSI was varied until all angular, energy and missing 
mass distributions of the input to the simulations where 
consistent with the data (see e.g. fig. |H1). From the com- 
parison of different simulations and the good agreement 
between the analysis of the three-photon events from the 
first measurement and the four-photon events from the 
second measurement of double 7r° production (see fig.l2U|l 
we estimate a systematic uncertainty of 5%. 
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Fig. 8. Missing mass spectra for double tt*^ photoproduction 
from Ca. Left side: incident photon energy of 500 MeV (be- 
low 77 threshold), Right side: incident photon energy 700 MeV 
(above 77 threshold). Closed circles: experiment, dashed his- 
tograms simulation of 27r° photoproduction, dotted histogram: 
simulated contribution of ri photoproduction, solid histogram: 
sum of both. 



4 The BUU-model 

In this section we discuss the main ingredients of the BUU 
transport model. For the present calculations we use the 
version discussed in detail in 34,35 . 

The model is based on the BUU equation which de- 
scribes the space-time evolution (r: space coordinate, p: 
momentum) of the spectral phase-space density Fi of an 
ensemble of interacting particles of type i = N, ^33(1232), 
TT, 77,... with mass /i: 

d \ 

— + VpH■Vr-VrH■VpjF^ir,p,^l;t)^ 



space density into an ensemble of discrete test particles 
propagating according to Hamilton's equations of motion. 

Before the actual reaction, the nucleon test particles 
are initialized in the nucleus according to a Woods-Saxon 
density function in coordinate space. In momentum space, 
we make use of the local density approximation, distribut- 
ing the nuclcons homogeneously in the Fermi sphere corre- 
sponding to the local density. For the baryons, the mean- 
field potential S is determined as described in 35 . In this 
work, we focus on a momentum dependent potential cor- 
responding to an equation of state with medium compress- 
ibility. 

The photon-nucleus reaction is modeled in terms of 
the absorption of the photon on a single nucleon, leading 
to final states P33(1232),Di3(1520),Sii(1535),_Fi5(1680), 
Nn, Ntttt, NV {V=p, w, 0), KA, KS and KKN. These 
entry states are selected with probabilities corresponding 
to the different reaction cross sections. They are then sub- 
jected to the final state interactions, which are described 
by the set of BUU equations. The cross sections for the 
7A reaction are then determined by averaging over an en- 
semble of such elementary reactions as outlined in |20| . 

Besides Fermi motion, binding effects and Pauli block- 
ing we also account for coUisional broadening of the most 
important resonances P33( 1232), Di3(1535) and Sii(1535). 
For the in-medium width of the P33 we either use the re- 
sults of Oset et al. 36 or the phenomenological spread- 
ing potential of Hirata et al. |37|, both giving roughly 
rcoii = 80 MeV at saturation density po- The latter is in- 
terpreted as fully absorptive, whereas the absorption part 
in 1221 is only about 40 MeV at po- For the Di3(152G) 
and the Sii(1535), we make use of the results of the self- 
consistent resonance-hole model in 0. In particular, for 
the D13 a total collision width of 2GGp//9o MeV at density 
p is used, while the Sn collision width is rather small and 
amounts only to 30p/ po MeV. In both cases, one quarter 
of the total width is used for the absorptive channels. 



Icoii[Fn,F^,Fp^3(i232),Fjj,...] (6) 5 Results and discussion 



The left-hand side - the Vlasov term - describes the prop- 
agation of the particles under the influence of a Hamilton 
function H , given by the expression 



H = ^i^, + sr + p^', 



(7) 



which in the case of baryons contains an effective scalar 
potential S [3S] . The right-hand side of the BUU equation 
(called the collision integral) consists of a gain and a loss 
term for the phase space density Fi at the different space- 
time points, accounting for interactions between the parti- 
cles beyond the mean-field potential. The collision integral 
contains collision rates for the different reaction types such 
as baryon-baryon and baryon-meson collisions, resonance 
formation and decay. They include cross sections for these 
processes and Pauli blocking factors for outgoing fermions. 
Due to the reactions of different particle types we end up 
with a set of coupled BUU equations, which is solved by 
a test particle ansatz. This converts the continuous phase 



5.1 Inclusive Tr'^-photoproduction 

The total inclusive tt^ cross section <7inci for the heavy nu- 
clei is compared in fig.^to the deuteron cross section. The 
cross sections are scaled by A^ff which throughout the pa- 
per is defined as Aeff=2 for the deuteron and A^jf = A^l'^ 
for the other nuclei. In this way, the average nucleon cross 
section is compared to the nuclear cross sections scaled 
in proportion to the nuclear surfaces. The shape of the 
total cross sections for the heavy nuclei is similar, but it 
differs from the deuteron. The peak in the A resonance 
region is broader for the heavier nuclei and shifted to a 
higher energy. Furthermore, the valley between the first 
and second resonance region is much more pronounced for 
the deuteron. The different shape is not explained by the 
larger Fermi momenta of the nucleons bound in the heav- 
ier nuclei. This is demonstrated with the dashed curve in 
fig. El which represents the deuteron cross section folded 
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with the nucleon momentum distribution of Ca. The latter 
was parameterized as ^^: 



PF{p)dp = K-p^ ■ (e"P'/" + ce-P/'') dp 



(8) 



with a = 0.42 fm^, 6 = 0.23 fm, and c=0.04. The folding 
was done as discussed in 03- Here, it was neglected that 
the deuteron excitation function already incorporates the 
(much smaller) Fermi smearing effects from the deuteron 
wave function which in principle would have to be de- 
folded before folding with the Ca distribution. 

It is tempting to argue that the shift of the peak in the 
A region by «20 MeV is related to an in-medium modifi- 
cation of the A resonance due to the A - nucleus interac- 
tion. Analyses of coherent tt" photoproduction gave indi- 
cations of a Z\ self-energy corresponding to a similar peak 
shift, almost independent of the nuclear mass 39,40,31 . 
However, the interpretation is not so simple. The shape 
of the low energy side of the peak reflects the superposi- 
tion of contributions from coherent and quasifrec single tt*^ 
photoproduction which have different energy dependen- 
cies. Furthermore, as discussed below, contributions from 
multi-body absorption effects of the incident photons are 
important in the high energy tail of the A. Altogether, 
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Fig. 9. Scaled total cross section aind of inclusive vr pho- 
toproduction from deuterium, carbon, calcium, niobium, and 
lead versus incident photon energy (scale factor 1/2 for the 
deuteron, A^'^ for all other nuclei). The solid line is only to 
guide the eye, the dashed line corresponds to the deuteron data 
folded with the momentum distribution of nucleons in calcium. 
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Fig. 10. Comparison of the inclusive cross section Oinci to 
BUU model calculations for calcium and lead. Open symbols: 
inclusive reaction, filled symbols: difference of inclusive and 
coherent reaction. Curves: BUU without FSI (dotted), BUU 
with FSI, A in-medium width from '36' (dashed) or |37) (solid). 



position and shape of the A peak are influenced by sev- 
eral effects of different nature and not simply connected 
to the in- medium properties of the A. 

The data for calcium and lead are compared to the 
BUU calculations in fig. EH Since by construction the 
BUU model accounts only for the incoherent part of pho- 
toproduction, the data are also shown after subtraction of 
the contribution from coherent single 7r° production taken 
from IHTI. The figure demonstrates the strong absorptive 
character of the FSI. The calculation without FSI overes- 
timates the lead data almost by a factor of four, but the 
inclusion of FSI produces a reasonable description of the 
data. The two models for the P33 in-medium widths lead 
to similar results. The inclusion of the spreading potential 
gives a larger suppression of the cross section. The broad 
structure of the first resonance peak cannot be reproduced 
by the values used for the P33 collision width but an in- 
crease much beyond 100 MeV at density po would lead to 
an even stronger suppression of the cross section and is not 
justified by theoretical models. However, it is known 021 
that two-body absorption mechanisms like 7NN— > NA are 
non-negligible in this energy range. Such contributions are 
not included in the BUU model. 

Typical angular distributions are summarized in fig. 
111! The contribution from coherent single tt" photopro- 
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duction is clearly visible for incident photon energies up 
to the A resonance region. It disappears at higher photon 
energies where the distributions are similar for all nuclei. 
The forward peaking at high photon energies is a conse- 
quence of the frame choice (photon - nucleus cm system). 
Due to the lack of the coherent part, the model calcula- 
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Fig. 11. Angular distributions in the photon - nucleus cm 
frame scaled by the mass number. Full symbols; inclusive cross 
section this work, open symbols; coherent tt" production 31 . 



tions cannot reproduce the low energy angular distribu- 
tions. Subtraction of the coherent part would not help, 
since as discussed for the deuteron in [T^IHT] the coherent 
and quasifree parts are somewhat interconnected via nu- 
cleon final state interaction. The effect is also observed for 
the heavy nuclei. As an example the angular distributions 
of the quasifree part around incident photon energies of 
330 MeV are shown in fig. ^1 The angular distributions 
are significantly depleted at forward angles where some 
strength is re-distributed to the coherent part. 

The momentum distributions of the pions are sum- 
marized in fig. 1131 These spectra show the influence of 
FSI effects. Apart from two body effects in the absorp- 
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Fig. 12. Angular distributions without coherent 7r° photopro- 
duction in the photon - nucleus cm frame scaled by the mass 
number for incident photon energies from 325 - 335 MeV. 



tion of the photons, the distributions from the deuteron 
can be regarded as source spectra of the pions, which are 
then modified for the heavier nuclei by FSI. At the lowest 
incident photon energies spectra from the deuteron and 
the heavy nuclei are similar. Up to energies of 500 MeV 
the distributions measured on the deuteron simply shift to 
larger momenta corresponding to the quasifree kinematics 
of single tt'^ production. However, the nuclei in this energy 
region show already a re-distribution of the pion momenta 
to smaller values. At small momenta the nuclear data scale 
with the mass number A which indicates a large mean 
free path of the pions in this range. At even higher inci- 
dent photon energies the pion spectra from the deuteron 
develop a double bump structure. The additional peak at 
lower momenta comes from double pion production reac- 
tions. At the highest energies this contribution becomes 
dominant. The contribution from single pion production 
for the heavier targets is suppressed due to FSI effects 
even at lower incident photon energies. The data for cal- 
cium and lead are compared to the BUU results obtained 
with the spreading potential |37] in fig. E| The discrep- 
ancy at the lowest incident photon energies is mainly due 
to coherent tt" photoproduction which contributes to the 
data but is not included in the BUU calculations. The 
general trend of the FSI effects, i.e. the shift of strength 
to small pion momenta, is reproduced by the model. 

The FSI effects are closely related to the absorption 
properties of nuclear matter for pions, that is the mean 
free path of pions in nuclear matter. The momentum de- 
pendence of this effect is analyzed in fig. El The main plot 
shows the ratio of the cross sections from the three heav- 
iest nuclei with respect to carbon, normalized to A^^"^. It 
is defined via: 



R, 



C2/3 



a(i2C)/122/3 



(9) 



The insert shows the scaling exponent a obtained from 
fits of the mass dependence of the data with the ansatz: 



a{A) oc A" 



(10) 



The differential cross sections are averaged over the entire 
range of incident photon energies from 200 - 800 MeV. The 
scaling coefficient is also shown for incident photon en- 
ergies above 400 MeV where contributions from coherent 
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Fig. 13. Distributions of the pion lab momenta for different in- 
tervals of incident photon energy. Left hand side: heavy nuclei, 
right hand side: deuterium. 



tt" photoproduction are negligible. Coherent contributions 
could obscure the scaling, since their production (before 
FSI) is not simply proportional to the mass number A. 
However, the plot shows that their influence is small. Both 
figures demonstrate qualitatively the same behavior. For 
pion momenta larger than roughly 200 MeV/c the cross 
sections scale like A"^^^ which indicates strong absorption 
and small mean free paths of the pions. At smaller pion 
momenta the nuclei become more transparent for pions 
and the cross sections scale with the volume. This is ex- 
pected since pions with a momentum of 227 MeV/c can 
excite a nucleon to the nominal mass of the A. Therefore, 
at momenta larger than Ri200 MeV/c pions are very likely 
re-absorbed via ttA^ — > A. This explains also qualitatively 
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Fig. 14. Distributions of the pion lab momenta for calcium 
(left hand side) and lead (right hand side) compared to BUU 
calculations. Data and model results are averaged over 100 
MeV wide bins of incident photon energies, the central val- 
ues are given in the plot. The calculation used the spreading 
potential of Hirata et al. 



the shift of strength to small pion momenta in fig. 1131 Only 
those pions which arc produced with small momenta or ac- 
quire small momenta after absorption and rc-cmission pro- 
cesses are likely to escape from the nucleus. Contributions 
from larger pion momenta stem mainly from the nuclear 
surface region. These momentum distributions are there- 
fore very sensitive to the FSI effects and thus a valuable 
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Fig. 15. Scaling of the momentum dependence of the differen- 
tial cross section with the nuclear mass number A. The main 
plot shows the ratio Rc2/3 (see eq. |3J for Ca, Nb and Pb. 
The insert shows the fitted coefficient a (see eg. llUH averaged 
over two ranges of incident photon energy. Curves: BUU results 
(solid: spreading potential |37| . dashed: Oset parameterization 
|3ti|). Both averaged over 400 - 800 MeV photon energy. 



testing ground for model predictions. The BUU calcula- 
tions show qualitatively similar behavior, however they do 
not reproduce the strong peak at pion momenta around 
100 MeV/c. This could indicate that the mean free path 
for low pion momentum is underestimated. However, one 
should be aware that the coefficient a does not only de- 
pend on the mean free path of the pions. It can be also 
influenced to some extent by the re-distribution of the pion 
momenta due to the absorption and re-emission processes. 
We have already noted in the discussion of the total 
inclusive cross section that in case of the heavier nuclei 
the valley between the first and second resonance region 
is much shallower than for the deuteron. This could be 
due either to FSI effects or to the initial photo-pion pro- 
duction mechanism. The two possibilities can be tested by 
a comparison of the momentum spectra to the deuteron 
data. For this purpose the ratio: 



i?rf = 






(11) 



which compares the nuclear data to the average nucleon 
cross section is plotted in fig. ll6l as function of the pion mo- 
menta for different regions of incident photon energies. At 
a given pion momentum the absorption probability should 
be independent of the incident photon energy. For the low- 
est incident photon energies the ratio shows the expected 
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Fig. 16. Ratio of cross section from heavy nuclei and deu- 
terium. The ratio Rd (see eq. 111^ for carbon, calcium, niobium 
and lead is plotted. 



behavior which reflects the momentum dependence of the 
pion mean free path discussed above: It is large for small 
pion momenta and then drops to unity around p^^ ~ 200 
MeV/c. However, at incident photon energies between 500 
- 600 MeV the ratio is always larger than unity, in partic- 
ular for momenta in the range between 200 to 350 MeV/c. 
This is an indication that for the heavy nuclei in the in- 
termediate energy range a certain fraction of pions is in- 
deed produced via mechanisms which do not operate or 
are less important for the deuteron. Possible candidates 
are contributions from multi-body absorption mechanisms 
like jNN — > NA (see e.g. 021 )j in- medium effects of any 
of the involved partial reaction channels, or charge ex- 
change reactions which effectively transfer strengths from 
final states like tt+tt" to the inclusive 7r° channel. How- 
ever, the latter are a less likely cause since the BUU calcu- 
lations do not reproduce the additional strength at inter- 
mediate photon energies although they include the charge 
exchange reactions. 



5.2 Semi-exclusive reaction channels 

In the following we will decompose the inclusive cross sec- 
tion into partial reaction channels. The photoproduction 
of 77 mesons contributes to the inclusive tt" cross section 
via the 77 -^ Stt" and ij — > tt^tt+tt" decays. Since the rj is 
long lived, it always decays outside the nuclei, so that these 
pions do not undergo FSI effects. The contribution to pion 
production is obtained directly from the measured 27 de- 
cay channel of the 77 mesons and the respective branch- 



ing ratios. The scaled total cross sections ari for rj pho- 
toproduction from nuclei are compared in fig. 1171 to the 
deuteron data |43| . In case of carbon and niobium previ- 
ous data from JH| are used. For calcium and lead results 
with much improved statistical quality, obtained with the 
improved setup, are shown. The nuclear data follow an 
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average cross section from the nucleon {ad/2) |43|. 



A^/^ scaling and agree above the production threshold on 
the free nucleon (Ethr=706 MeV) with the average nu- 
cleon cross section. In the immediate threshold region the 
reaction is only possible for large momenta of the bound 
nucleons anti-parallel to the incident photon momentum. 
The reaction on the deuteron is therefore suppressed in 
the threshold region since the Fermi momenta are much 
smaller than for the heavier nuclei. 

Quasifree single tt*^ photoproduction was extracted with 
the kinematical cuts discussed in sec. 13.21 Results for the 
second resonance region were already published in |24| . 
Again the nuclear cross sections a'^l scale like A^^^ (see fig. 
I18|) and are similar to the average nucleon cross section. 
The difference between the scaled nuclear and deuteron 
cross sections at low incident photon energies is related to 
coherent tt" photoproduction which is more important for 
the heavier nuclei. Note that in this channel in contrast to 
the inclusive reaction the position of the A peak is shifted 
to a lower incident photon energy for the heavier nuclei. 

In fig. El the data for calcium and lead arc compared 
to the results of the BUU model. Again, the data are also 
plotted after subtraction of contributions from coherent 
tt'^ production to facilitate a better comparison with the 
model. The calculations include the same missing energy 
cuts as the experiment. The curves show a similar behavior 
as in the case of the inclusive cross section in fig. El In par- 
ticular, the discrepancy at intermediate photon energies 
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Fig. 18. Total cross section a'^l of single quasifree tt" photo- 
production from heavy nuclei. The nuclear cross sections are 
normalized to A^ " and compared to the average cross section 
from the nucleon {ad/2). The insert compares for ''"Ca a'^l (full 
points) to the cross section with subtracted coherent part. The 
curve corresponds to the deuteron cross section {ad/2) folded 
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also shows up here. This suggests that two-body photon 
absorption mechanisms also contribute to this channel. 
In the second resonance region, the calculations overesti- 
mate the data. However, this does not necessarily mean 
that this is solely due to an overestimate of the D13 contri- 
bution which in the calculations only amounts to roughly 
25 % of the total cross section for E-, = 700 MeV. 

The experimental results for double tt*' photoproduc- 
tion are summarized in fig. 1201 {a2T,o ) . In the first mea- 
surement the limited solid angle did not allow an analysis 
of events with four photons. Those data are therefore lim- 
ited to the practically background free energy range below 
the 77-production threshold (triangles in the figure for C, 
Ca, Nb, Pb). The measurement with the improved setup 
allowed the clean extraction of the cross section from 4- 
photon events via a missing mass analysis. Background 
from 77 decays was subtracted. The results for calcium and 
lead (filled circles) are shown in the figure up to a max- 
imum photon energy of 800 MeV. All nuclear data scale 
like A^/^ and agree with the average nucleon cross section 
crd/2 taken from |;-{2| . 

The same result is found for tt'^tt* production which is 
shown in fig. [21 The data for calcium and lead again scale 
like A^/'^. Also shown in this figure are the cross sections 
for jp — > mr'^TT'^ 0^ and 777. -^ ptt^tt" 0^1 ■ Finally, the 
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cross section for 70? — > npir^TT was obtained from the 
inclusive tt" photoproduction cross section Ginci on the 
deuteron by subtraction of the cross sections from all other 
partial channels (see eq. ^ . It is obvious from the figure 
that also in this case the nuclear cross sections scaled by 
^^2/3 agj-gg with the average nucleon cross section. 

The results for the semi-exclusive quasifree reaction 
channels thus can be summarized in the following way: all 
investigated nuclear cross sections are related to a good 
approximation to the deuteron cross section by: 






(12) 



The scaling among the heavier nuclei holds even when 
the comparison to the deuteron is disturbed, e.g. by the 
effects of coherent single 7r° production at low photon en- 
ergies and or by effects of Fermi smearing close to the 77 
threshold. The A^^^ scaling is the limiting case of strong 
FSI effects due to the short pion mean free path. This 
means that the quasifree exclusive reactions probe only 
the nuclear surface region. Since the properly scaled nu- 
clear cross sections agree with the deuteron cross section 
no significant in-medium effects are observed for the low 
density surface zone of the nuclei. 
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Fig. 20. Total cross section 0-2^0 of double tt" photoproduction 
normalized to A^'^ and compared to the average nucleon cross 
section {od/2). 



With the knowledge of the nuclear cross sections for 
7] and double 7r° photoproduction it is possible to correct 
the contribution from higher pion multiplicity channels 
to the total inclusive cross section and thus to determine 
the part of the total photoabsorption cross section (cr„m , 
see eq. EJ which corresponds to events with at least one 
TT^ or one rj in the final state. The result is compared 
in fig. 1221 to the BUU model calculations. Again also the 
results with subtracted coherent part are shown. The cal- 
culations show essentially the same behavior as already 
observed for the inclusive cross section in fig. 9. In par- 
ticular for the lighter nuclei the second resonance peak 
is still overestimated and the valley between the A reso- 
nance and the second resonance region is too pronounced. 
In the case of carbon it should be noted that the assump- 
tions of the BUU model are much less fulfilled for such 
a light nucleus with a relatively large surface/volume ra- 
tio. However, a further observation is that while for the 
heaviest nuclei the agreement in the second resonance re- 
gion is much improved, the disagreement in the 'valley 
region' around 500 MeV grows. This could mean that the 
agreement in the second resonance region is only acciden- 
tal because the contributions from two-body absorption 
effects are neglected. A quantitative investigation of the 
influence of the two-body absorption effects, in particular 
the dependence on the incident photon energy, is therefore 
highly desirable. 
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6 Summary and conclusion 

Cross section data for meson photoproduction from nuclei 
in the energy range from 200 - 800 MeV have been ob- 
tained for all final states which involve at least one neutral 
meson. The inclusive data have been used for an analysis 
of the pion final state interaction in nuclei. The momen- 
tum dependent absorption properties of nuclear matter 
show qualitatively the expected behavior: strong absorp- 
tion for pions which can excite the A resonance on the 
nucleon and almost perfect transparency of the nuclei for 
low energy pions. The absorption properties within the 
BUU model show a similar behavior, although no quanti- 
tative agreement is achieved for small pion momenta. 

The cross sections of all exclusive quasifree meson pro- 
duction reactions scale with the nuclear surface {A^/"^) and 
agree to a good approximation with the respective cross 
sections for the deuteron normalized to the deuteron mass 
number. This means that no significant in-medium effects 
are observed for the quasifree processes which test the 
low density surface region of the nuclei. Even the inclu- 
sive data, which include some fraction of non-quasifree 
reactions, show a clear signal for the second resonance 
bump. This is an indication that the complete suppres- 
sion of this bump observed in total photoabsorption reac- 
tions is not the result of trivial nuclear effects like Fermi 
motion. The observation that the excitation function for 
meson production from the surface region of nuclei has 
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Fig. 22. Comparison of the total cross section for neutral me- 
son production (cr„m, see eq. |2J to BUU model calculations. 
Open symbols: a„m, filled symbols: coherent single tt" produc- 
tion subtracted. Curves: BUU, A in-medium width from Oset 
parameterization (dashed) |36| or Hirata (solid) |37|. 



a different energy dependence than total photoabsorption 
in the nuclear volume seems to indicate density dependent 
in- medium effects. 

The experimental data have been compared in detail 
to calculations in the framework of the BUU model. It is 
observed that the model is remarkably succesful in repro- 
ducing FSI effects such as pion absorption and the shift of 
strength in the pion momentum spectra. However, some 
discrepancies on the quantitative level remain. In the A- 
resonance region the inclusion of realistic P33 in-medium 
widths lead to a better description of the data. Also in the 
second resonance region the calculations are close to the 
data although especially for lighter nuclei an overestima- 
tion is visible. The intermediate energy region is not repro- 
duced which is in line with the theoretical findings in pho- 
toabsorption that multi-body photon absorption mecha- 
nisms, which are not included in the BUU model, play an 
important role in this region. Effenberger et al. 46 have 
calculated such effects for photoabsorption on '^''Ca. They 
find them negligible below incident photon energies of 300 
MeV and rising to a 25 % effect at photon energies around 
500 MeV. A comparable contribution to pion production 
would much improve the agreement between data and cal- 
culations in this energy region. However, at present it is 
not known how much the two-body effects will contribute 
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at higher incident photon energies. A quantitative treat- 
ment of these effects is thus necessary in future theoretical 
studies of the in-medium effects in the second resonance 
region. 
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